Double deeply virtual Compton scattering off the nucleon 
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We study the double deeply virtual Compton scattering (DDVCS) process oflt the nucleon, through 
the scattering of a spacelike virtual photon with large virtuality resulting in the production of a 
timelike virtual photon, decaying into an e"'"e~ pair. This process is expressed in the Bjorken regime 
in terms of generalized parton distributions (GPDs) and it is shown that by varying the invariant 
mass of the lepton pair, one can directly extract the GPDs from the observables. We give predictions 
for the DDVCS cross section and beam helicity asymmetry and discuss its experimental feasibility. 

PACS numbers: 13.60.Fz, 12.38.Bx, 13.60.Le 



The understanding of hadron structure in terms of 
quark and gluon degrees of freedom, remains an out- 
standing challenge. An important source of information 
is provided by experiments involving electroweak probes. 
In this way, elastic form factors as well as quark and 
gluon distributions in the nucleon have been mapped out 
in quite some detail. In recent years, a whole new class 
of hard exclusive reactions have become accessible both 
theoretically and experimentally to study hadron struc- 
ture. In particular, the deeply virtual Compton scatter- 
ing (DVCS) and hard electroproduction of meson pro- 
cesses are at present under investigation at different fa- 
cilities (HERMES §, JLab §, HERA |, |), or wiU be 
addressed by experiments in the near future. In these 
processes, a highly virtual photon (with large virtual- 
ity Q^) scatters from the nucleon and a real photon (in 
the case of DVCS) or a meson is produced. Due to the 
large scale involved, these hard exclusive processes 
are factorizable in a hard part, which can be calculated 
from perturbative QCD, and a soft part, which contains 
the information on nucleon structure and is parametrized 
in terms of generalized parton distributions (GPDs) (see 
Refs. [|[ H for reviews and references therein). 

The GPDs depend upon the different longitudinal mo- 
mentum fractions x + ^ (x — ^) of the initial (final) quarks 
(see upper left panel of Fig. As the momentum frac- 
tions of the initial and final quarks are different, in con- 
trast to the forward parton distributions, one accesses in 
this way quark momentum correlations in the nucleon, 
which are at present largely unknown. Furthermore, sum 
rule integrals of GPDs over x provide new nucleon struc- 
ture information and are also amenable to lattice QCD 
calculations for direct comparison. In particular, the sec- 
ond moment of a particular combination of GPDs gives 
access to the total angular momentum carried by quarks 
in the nucleon . Such a quantity would be highly com- 
plementary to the information extracted from polarized 
deep-inelastic scattering experiments, which found that 
about 20 - 30 % of the nucleon spin originates from the 
quark intrinsic spins (see Ref. |9j for a recent review). 

To obtain these new informations, one of the main chal- 
lenges is to directly extract the GPDs from observables. 



In the DVCS or hard exclusive meson electroproduction 
observables, the GPDs enter in general in convolution in- 
tegrals over the average quark momentum fraction so 
that only ^ (half the difference of both quark momentum 
fractions) can be accessed experimentally. A particular 
exception is when one measures an observable propor- 
tional to the imaginary part of the amplitude, such as 
the beam helicity asymmetry in DVCS. Then, one actu- 
ally measures directly the GPDs at some specific point, 
X — ^, which is certainly an important gain of informa- 
tion but clearly not sufhcient to map out the GPDs inde- 
pendently in both quark momentum fractions, which is 
needed to construct sum rules. In absence of any model- 
independent "deconvolution" procedure at this moment, 
existing analyses of DVCS experiments have to rely on 
some global model fitting procedure. 






FIG. 1: Diagrams for the Ip lpe'^e~ process : DDVCS 
process (upper left), vector meson (VM) production process 
(upper right), Bethe-Heitler (BH) processes (lower two dia- 
grams). Crossed diagrams are not shown but also included. 



The double DVCS (DDVCS) process, i.e. the scatter- 
ing of a spacelike virtual photon from the nucleon with 
the production of a virtual photon in the final state, pro- 
vides a way around this problem of principle. Compared 
to the DVCS process with a real photon in the final state, 
the virtuality of the final photon in DDVCS yields an 
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additional lever arm, which allows to vary both quark 
momenta x and ^ independently. This additional lever 
arm of the DDVCS compared to the DVCS process has 
already been noted in Refs. |10, 1^, Also, the lepton 
pair production process induced by a real photon, i.e. the 
7P e^e^p reaction has been studied 1 13|. In this letter, 
we provide the first numerical estimates of the DDVCS 
and its competing processes and show how the GPDs can 
be directly extracted from DDVCS observables. 

The DDVCS process can be accessed through the 
Ip —> lpe^e~ reaction, which is characterized by the four- 
momenta : k [k') of the incoming (scattered) leptons /, 
p {p') of the initial (final) nucleons, and /~, l'^ of the lep- 
tons in the produced e+e~ pair. 

To see how the DDVCS process can yield more com- 
plete information of GPDs than the DVCS process, one 
first has to discuss its richer kinematics. The DDVCS 
process is characterized by 8 independent kinematical 
variables. Firstly, there are the same 5 kinematical vari- 
ables which specify the DVCS process and which we 
choose as : the initial beam energy Eg] the virtuality 
of the incoming photon in the upper left diagram of 
Fig. ^ i.e. = —q"^, where q = k — k'; the usual Bjorken 
variable xb = / (2p.<z); the four- momentum transfer to 
the nucleon t = A^, where A = p' — p; and the out-of- 
plane angle $ between the production plane, spanned by 
the vectors q and q' , and the scattering plane spanned 
by the vectors k and k'. Furthermore, one needs 3 addi- 
tional variables to fully characterize the DDVCS process 
which we choose as : the virtuality q'^ = {l^- + l^+Y of 
the produced e^e~ pair; and the 2 angles of one lepton of 
the produced lepton pair, evaluated in the cm. system 
of the e^e~ pair, and which span the solid angle di^*- . 

At large Q^, we calculate the DDVCS process in the 
handbag approximation as shown in Fig. ^ (upper left 
diagram), which yields the amplitude : 
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where /i {v) refer to the four-vector indices of the incom- 
ing spacelike (outgoing timelike) virtual photons respec- 
tively, n is a light-like vector along the direction of the 
incoming virtual photon, and where we r efer to Ref. § 
for the expressions of the symmetrical (antisymmetri- 
cal) twist-2 tensors {e^^ ). Furthermore, in Eq. (|^), 



N{p), N{p') represent the nucleon spinors and tun is the 
nucleon mass. The GPDs H, E, H, E in Eq. (^ are the 
same as in the DVCS case, and depend on the arguments 
x, ^, and i, with x and ^ as defined in Fig. |^. The coef- 
ficient functions in the DDVCS amplitude of Eq. (|l|) 
take the form : 



1 



± 



1 



(2^' + x + (2e - 
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where —2^' and 2(f — ^') are the longitudinal momentum 
fractions of the incoming spacelike and outgoing timelike 
virtual photons respectively (see Fig. In the large 
Q2 limit, one has 2^' — > xb/{^ — xb/'^)- The difference 
(2^' — ^) appearing in the quark propagators in Eq. (|^) 
can be expressed as (relative to ^) : 



2^'-e 



l-((7'2-A2)/g2 + 8e'2mVg2 



1 + (g'2 - A2)/Q2 



1 - q'VQ^ 
1 + g'VQ' 



(3) 



with — 



- A2/4. For the DDVCS process, by vary- 



ing the virtualities of both incoming and outgoing virtual 
photons, one can vary independently the variables, ^ and 
whereas, in DVCS, only one variable can be varied as 
S, ^ C ■ One then sees from Eqs. (|l|,||) that the imaginary 
part of the DDVCS amplitude (which can be directly 
measured through the beam helicity asymmetry as dis- 
cussed further on) will access, in a concise notation, the 
GPD(2^' - C,C,i), and allows to map out the GPDs as 
function of its three arguments independently. In the 
second line in Eq. (^, we have indicated the expression 
in the large limit, which is displayed in Fig. 0. For 
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FIG. 2: Range in the argument x — 2^'— ^ (relative to ^) of the 
GPD {x, ^, t) which one accesses by measuring the imaginary 
part of the DDVCS process at different values of the lepton 
invariant mass q'^ (relative to the initial virtuality Q^). 

a timelike virtual photon (i.e. q'^ > 0), one can only 
access the a; < ^ region in the arguments {x,^,t) of the 
GPDs because, for kinematical reasons, q'^/Q^ will al- 
ways be less than 1. Therefore, the imaginary part of 
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the DDVCS amplitude maps out the GPD where its first 
argument Hes in the range < 2^' — ^ < ^. In particular, 
when q'"^ is varied from to (5^/2, the argument x spans 
about 2/3 of the range [x,^]. Although one does not 
access the whole range in x, clearly, the gain of informa- 
tion on the GPDs is tremendous as no deconvolution is 
involved to access this region of the GPDs. Furthermore, 
a; < ^ is just the range where the GPDs contain wholly 
new information on mesonic {qq) components of the nu- 
cleon, which is absent in the forward limit (where ^ = 0). 
To access the range a; > ^ one would need two spacelike 
virtual photons, necessitating to select the two-photon 
exchange process in elastic electron nucleon scattering. 

Besides the DDVCS process, the Ip lpe^e~ reaction 
contains two classes of Bethe-Heitler (BH) processes as 
shown in Fig. ^ The BH processes are fully calculable as 
they involve elastic nucleon form factors. Furthermore, 
the outgoing timelike photon which couples pointlike to 
the quark line in the DDVCS process can also originate 
from a neutral vector meson (VM) which couples to the 
quark line through a onc-giuon exchange (upper right di- 
agram in Fig. . For the contamination of the VM pro- 
duction, we estimate it by the leading order amplitude 
for the hard electroproduction of longitudinally polar- 
ized VM iQ. For this process, which is of order 0{as) 
in the strong coupling constant compared to the handbag 
process, a factorization theorem has been proved al- 
lowing to express its amplitude also in terms of GPDs. 

In the following, we will estimate the coherent sum of 
all these processes. The fully differential cross section of 
the Ip lpe^e~ reaction can then be expressed as : 
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where y = {p.q)/{p.k), and where Tbh, Tddvcs, and 
TvM are the amplitudes for the BH, DDVCS and VM 
processes respectively. When integrating Eq. (Q) over the 
angles of the produced e+e^ pair, the resulting DDVCS 
cross section reduces in the limit q'^ ^ to : 



da 



dQ"^ dxB dtd^dq''- 
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dQ^ dxsdtd^ 



N 

772 



(5) 



where the DVCS cross section appears on the rhs of 
Eq. (^. The factor N in Eq. (||) is given by 



N = 



An 



A 

3' 



(6) 



where ae ~ 1/137 is the fine structure constant, intro- 
duced by the decay of the outgoing photon into the lepton 
pair. One sees that the downside of the DDVCS process 
is that it involves small cross sections : at a virtuality 



q'^ = 1 GeV^, the DDVCS cross section is reduced by at 
least a factor 7V"^ (w 1.3 x 10^) compared to the DVCS 
cross section. At lower values of q''^, the DDVCS cross 
section rises however as l/q'^- 

Besides the DDVCS cross section, a particularly infor- 
mative observable is obtained by scattering a longitudi- 
nally polarized lepton beam and flipping its helicity. The 
resulting single spin asymmetry (SSA) originates from 
the interference of the DDVCS and BH processes as : 



SSAr^lm [Tbh {Tddvcs + Tvm)*] 



(7) 



Because the BH process is real, the SSA accesses the 
imaginary part of the DDVCS -I- VM process, which is 
proportional to the GPD(2^' - ^,^,t) (see Eqs. (|J2|)). 
In Fig. ^, we show the q'^ dependence of the estimated 
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FIG. 3: Cross section (upper panel) and SSA (lower panel) of 
the ep — > epe^e~ process as function of the e^e~ virtuality 
q'^. Dashed curves : BH processes; dashed-dotted curves : 
BH + DDVCS processes, full curves : BH + DDVCS + pi 
processes. The dotted curves are the corresponding results for 
the ep — > ep7 proces. The ep — > epe^e" cross section is scaled 



with iV" 



q , in order to reproduce exactly the e; 



cross section in the limit q'^ 0, according to Eq. ( 



ep7 



cross section and SSA for the ep — )■ epe~^e~ process in 
kinematics accessible at JLab. As the twist-2 SSA basi- 
cally displays a sin $ structure, we show its value at $ = 
90°. For the GPDs, we use a ^-dependent parametriza- 
tion (see Refs. 0, |4))' ^^i^S MRSTOl |l§ forward 
quark distributions as input. As is seen from Fig. |^, 
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epe+e 



^g'^ reduces to the 



we firstly confirm numerically that the ep 
cross section scaled with the factor TV 
ep — > ep7 cross section when approaching the real photon 
point. Similarly, the SSA for the ep — > epe^e~ process 
reduces to the corresponding SSA for the ep epj pro- 
cess. When going to larger virtualities q'^, the DDVCS 
shows a growing deviation from the 1/q'^ behavior and 
the magnitude of the SSA decreases. Furthermore, we 
show in Fig. |^ the contribution of the p^, ^ e+e^ pro- 
cess (upper right diagram of Fig. |l]), which is the most 
pronounced VM process. We find that, except in the im- 
mediate vicinity of q'^ « to^, the — > e^e" process is 
very small. This can be understood because the cross 
section for the p^, ~^ e+e" process is reduced by a factor 
compared to the p° process, whereas the DDVCS 
cross section is only reduced by a factor compared 
to the DVCS cross section. Similarly, the SSA is only 
slightly affected by the VM process and is dominantly 
proportional to the imaginary part of the DDVCS pro- 
cess according to Eq. (j^) . The strong sensitivity of the 
SSA on q'^, as seen from Fig. ^, should therefore allow 
to map out the GPDs in the range a; < ^. 

Fig. ^ shows the comparison between the ep — > 
ep(7, ^iid ep —> ep(7, p^,) ~* ep{e^e~) processes for a 
typical kinematics accessible at JLab at 6 GeV. Whereas 
the ep —>■ epp'j^ process is roughly comparable to the 
DVCS -I- BH one for these kinematics, their "timelike" 
analogues, show that the p° channel is suppressed by 2 
orders of magnitude with respect to the DDVCS-f BH. 

Given that about lO'^ DVCS-^BH events were mea- 
sured recently at CLAS j^] in an effective 4-day data 
taking period at a luminosity of 10^^cm~^s~^ in a non- 
dedicated experiment, it can certainly be envisaged that 
the DDVCS-I-BH cross section, which is about 3 orders 
of magnitude lower, be measured with a dedicated long- 
time experiment. A luminosity of 10'^^cm~^s~^, pro- 
jected at CLAS for the upgrade of JLab at 12 GeV, or at 
a future dedicated lepton facility, would allow to measure 
this reaction with reasonable statistics. 

In conclusion, we have studied the DDVCS process 
with the production of a timelike virtual photon, decay- 
ing into an e+e^ pair. We have expressed the DDVCS 
amplitude in terms of GPDs and have shown that by 
varying the virtuality q'^ of the timelike photon, one can 
map out the GPDs as function of both initial and final 
quark momentum fractions. We have given cross section 
estimates for the DDVCS and its associated processes. 
Although the cross sections are small, their measurement 
seems feasible with a dedicated experiment at JLab and 
at a future high-energy, high-luminosity lepton facility. 
Of particular interest is the SSA using a polarized lepton 
beam. We have shown that by measuring the SSA, one 
can directly extract the GPDs in the domain where one 
is sensitive to qq correlations in the nucleon, providing a 
whole new source of nucleon structure information which 
is absent in forward quark distributions. 
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FIG. 4; Comparison of the different cross sections (as indi- 
cated on the curves) for the ep ep(7)PL) (left panel) and 
ep — » epe^e~ (right panel) reactions in JLab kinematics. 



This work was supported by the Deutsche Forschungs- 
gemeinschaft (SFB443), by the CNRS/IN2P3, and in 
part by the European Commission IHP program (con- 
tract HPRN-CT-2000-00130). 



[1 
[2: 

[s: 
[4: 

[5] 
[6] 

[7 

[8 

[9: 

[10: 

[11 
[12 

[13 
[14 

[15 
[16 



A. Airapetian et al. (HERMES Collaboration), Phys. 
Rev. Lett. 87, 182001 (2001). 

S. Stepanyan et al. (CLAS Collaboration), Phys. Rev. 
Lett. 87, 182002 (2001). 

C. Adloff et al. (HI Collaboration), Phys. Lett. B 517, 
47 (2001). 

F.R. SauU (ZEUS Collaboratio n). Proc. ICHE P 1999, 
(Tampere, Finland, July 1999), tiep-ex/0003030l 
X. Ji, J. Phys. G 24, 1181 (1998). 

A. V. Radyushkin, in the Boris loffe Festschrift 'At the 
Frontier of Particle Physics / Handbook of QCD', edited 
by M. Shifman (World Scientific, Singapore, 2001). 

K. Goeke, M.V. Polyakov, M. Vanderhaeghen, Prog. 

Part. Nucl. Phys. 47, 401 (2001). 

X. Ji, Phys. Rev. Lett. 78, 610 (1997). 

B. W. Filippone and X. Ji, Adv.Nucl.Phys. 26, 1 (2001). 

D. Miiller, D. Robaschik, B. Geyer, F.-M. Dittes, and J. 
Horejsi, Fortschr. Phys. 42, 101 (1994). 

J. Blumlein, D. Robaschik, Nucl.Phys. B 581 449 (2000). 
A.V. Belitsky, D. Miiller, and A. Kirchner, Nucl. Phys. 
B629, 323 (2002). 

E. R. Berger, M. Diehl, B. Pire, Eur. Phys. J. C 23, 675 
(2002). 

M. Vanderhaeghen, P.A.M. Guichon, and M. Guidal, 
Phys. Rev. Lett. 80, 5064 (1998); Phys. Rev. D 60, 
094017 (1999). 

J.C. Collins, L.L. Frankfurt, and M. Strikman, Phys. 
Rev. D 56, 2982 (1997). 

A.D. Martin, R.G. Roberts, W.J. Stirling, R.S. Thorne, 
Eur. Phys. J. C 23, 73 (2002). 



